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A discrete WKB method is developed for calculating tunnel splittings in spin prob- 
lems. The method is then applied to the issue of how nuclear spins affect the 
macroscopic quantum coherence of the total magnetic moment in small magnetic 
particles. The results are compared with numerical work, and with previous in- 
stanton based analytic approaches. 



1 Motivation for this Work 

Small magnetic particles have now been investigated as good candidates for 
experimental obserisation of macroscopic quantum tunneling (MQT) and co- 
herence (MQC).Erlra This is because at first sight the main critera for a system 
to be a good candidate for seeing MQP (P for phenomena) are met. One 
of these is that the energy barrier through which the system must tunnel be 
microscopic, even though the tunneling variable itself is macroscopic. This 
demand is met for the physical reason that the anisotropy energy barrier orig- 
inates in spin-orbit or spin-spin interactions at the microscopic level, both of 
which are relativistic effects, and hence small. A second criterion is that there 
be a well defined macrovariable, whose dynamics can, to a first approximation, 
be isolated from those of other microscopic degrees of freedom. This crite- 
rion can be met by working at temperatures sufficiently below the equivalent 
anisotropy energy gap, as spin wave excitations are then frozen out, and we 
may focus on the net magnetic moment of the particle as the macrovariable. 

There are also extremely good reasons to believe, however, that MQP 
are very hard to observe in general. Chief among these reasons is that the 
couplings of the macrovariable to the microscopic degrees of freedom give rise 
to decoherence. This effect can be especially severe in the case of macroscopic 
quantum coherence, even when the coupling is so small that the semiclassical 
dynamics of the macrovarible are significantly underdamped. The best studied 
example which cmvides a detailed illustration of this point is that of the spin- 
boson problem. qO'Q Further support for this point is provided by the scantiness 
of observational evidence for quantum coherence even in systems that one 
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would regard as microscopic. Thus, we know of only a few dozen or so "flexible" 
molecules like KH3 which display coherent flip-flop between different nuclear 
configurations. El The frequency of flip-flop is generally in the 0.1-100 GHz 
range. Given the vast range of molecular structures, and bonding strengths, 
there must surely be many naturally occurring molecules whose energy barriers 
and attempt frequencies are such as to put the flip-flop frerquency at anout 
1 Hz. Yet such flip-flop has never been seen. The reason almost certainly is 
that such molecules are never naturally encountered in isolation by themselves, 
and collisions and other environmental interactions are very effective in wiping 
out the quantum coherence. 

It is thus important to identify mechanisms for decoherence in the mag- 
netic particle system, and several have been put forth (phonons, magnons, 
Stoner excitations). The most critical, however, is the spin of the nuclei in the 
particle. The hyperfine coupling between the nuclear and electronic spins in 
magnetic solids is of order 100 MHz or more (in frequency units) per nucleus, 
which is rather high on the scale of the expected MQC frequencies. At the 
same time this frequency is rather low on the scale of the attempt frequencies 
associated with the electronic moments. Nuclear spins are therefore likely to 
be extremely efficient decoherers or "observers" of the direction of the moment 
of a small particle. This expectation is confirmed by theoretical calculations 
for both MQT,@ and MQC.u It is the latter that I wish to focus on and revisit 
in this article, for several reasons. First is that this d£caherence mechanism 
falls outside the scope of the harmonic oscillator bath,E9tj so one cannot rely 
on previous results. The calculations in Ref. 9 are done using an instanton 
technique with many physically motivated approximations about the trajecto- 
ries likely to give the dominant contribution to some path integral. It is not 
obvious even to me that this calculation is done in strict adherence with the 
Rheinheitsgebot. On the occasions that I have given seminars on the subject, 
the quizzical looks on the faces of my audience make it clear to me that it 
does not fully believe or understand my approach. It is therefore desirable to 
study this problem using different methods, and we shall do so in this article 
using the discrete WKB method. Although this method itself is quite old (see 
Braun's review i3 for references), it does not appear to have been usetLin spin 
tunneling problems except for some work by van Hemmen and Siito.Ej These 
authors have not fully exploited the power of this method, however. We shall 
see that calculations which have traditionally been done by instanton methods 
can be done much more simply to the same accuracy using this method. The 
opportunity to present the technical aspects of this method in a tutorial vol- 
ume devoted to tunneling in complex systems is greatly welcome, and provides 
me with another reason for writing this article. 
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It should be noted that the same subject has_also been studied by Prokof 'ev 
and Stamp in several papers. The first of theseliJ correctly notes that the tun- 
neling amplitude is suppressed by nuclear spins, but fails to recognize the 
implications of this fact for the nature pi. the tunneling spectrum. See footnote 
d for more on this. Subsequent papers cB puport to make detailed calculations 
of the tunneling spectrum including lineshapes, and also to include a host of 
other physical effects, such as spin diffusion, the Suhl-Nakamura interaction, 
and others with less familiar names such as "topological decoherence" , "or- 
thogonality blocking" , and "degeneracy blocking" . I cannot comment on the 
later work, simply because I do not understand much of it, especially some of 
the more mathematically specific and detailed conclusions, about lineshapes, 
for example. My goal in this article will be much more modest. It was argued 
in Refs. 3 and 9 that in the presence of nuclear spins the tunneling spectrum 
is broken into several resonance lines. A physical interpretation was attached 
to this broken-up spectrum, and formulas were presented for their frequencies 
and spectral weights. This article corroborates these claims. As in the pre- 
vious papers, I have not attempted to model or calculate the details of the 
relaxation, i.e., the lineshapes. To this extent, I am only prepared to claim 
a qualitative understanding for the very low frequency part of the tunneling 
spectrum. Fortunately, it is the high frequency end that is most likely to be 
experimentally relevant if at all, and here the situation is much better. 

The plan of this article is as follows. In Sec. 2, I will give a brief intro- 
duction to the physical problem and the models studied in this paper. Sec. 3 
contains a general discussion of the discrete WKB method. This is used to 
calculate the bare tunnel splitting, i.e., without nuclear spins, in Sec. 4. Nu- 
clear spins are added to the problem in Sec. 5. This section has the bulk of 
the new results in this article. The results for the tunnel splitting(s) obtained 
via the discrete WKB method are checked against those from exact numerical 
diagonalization of modcLHamiltonians. I also compare them with those from 
the instanton approach, □ and identify which features of the latter appear to 
be quantitatively robust, and which are only qualitatively correct. I conclude 
in Sec. 6 with a summary of the effects of nuclear spins, and some general 
remarks on the observability of MQC in magnetic particles. 

2 Introduction to Physical Problem and Models 

2.1 Physical System 

The physical system is a small insulating magnetic particle, about 50 A in di- 
ameter, at millikelvin temperatures. The anisotropy energy gap is much larger 
than this, so spin wave excitations are frozen out, and the individual atomic 
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spins are orientationally locked together. Furthermore, at this size the particle 
typically contains only one magnetic domain. Thus the magnitude of the total 
spin or magnetic moment of the particle is essentially fixed proportional to the 
number of atomic spins (assuming one magnetic species for simplicity), and 
the only relevant dynamical variable is its direction. (It may be useful to think 
of the system in terms of a Heisenberg-like model, with additional single-ion 
anisotropy terms.) The existence of anisotropy implies that some directions are 
energetically favored over others, and we wish to investigate whether the spin 
orientation can display quantum mechanical behavior, in particular tunneling 
and/or coherent oscillation between different energy minima. _ 
A model Hamiltonian which incorporates the above features is:EJ 

Ho = -KSl + k' 2 S 2 x (1) 

Here S x , S y , and S z are the components of our large spin, and fc- are phe- 
nomenological anisotropy coefheents. This particular Hamiltonian is time- 
reversal invariant, which should be the case if our particle is not subject to any 
external magnetic fields. We take k[ > 0, k' 2 > 0, so that in classical language, 
±z are easy directions, and ±x are hard directions. Quantum mechanically, 
the two classical ground states ±z will be split by tunneling. A generally valid 
approximate expression for the tunnel splitting can be obtained by writing 

A ~ uj e exp(-S , u/w e ). (2) 

[A more exact expression is given in Eq. (|24|) below.] In this formula, 

u e = 2{k 1 k 12 ) 1 ' 2 (3) 

is the oscillation (or precession) frequency for small deviations of the spin 
orientation from the classical equilibrium directions ±z,^ and we have also 
defined 

kx = Sk[, k 2 = Sk' 2 , k t 2 = ki + k 2 . (4) 

Further, Su is a quantity proportional to the energy barrier, of order Sk±2- 
We work throughout in units such that Ti = 1. 

We next wish to consider the influence of nuclear spins. To model this let 
us suppose that N n of the atoms have nuclei with spins Ij , and take all of these 
to be of magnitude 1/2. Let us further simplify the problem and assume that 
the hyperfine interaction for each of these atoms is of identical strength and of 
the form Sj ■ Ij, where Sj is the electronic spin on atom i. All the assumptions 

"This result may be derived by writing down the Heisenberg equation of motion for S, and 
linearizing it in small deviations from S = Sz. 
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except that of identical interaction strength are immaterial, and even this is a 
rather good approximation. The effects of relaxing it will be briefly discussed 
in Sec. 5.3. The total hyperfine contribution to the Hamiltonian can then be 
written as 

W hf = --Vs i -I 1 ; = -^S.I tot , (5) 

s b 
%—\ 

where we have expressed the coupling constant in terms of u> n , the nuclear 
Larmor frequency that would be obtained if the electronic spin orientation 
were fixed. The first expression in Eq. (|j) is just the sum of the interactions 
for the individual atoms, and the second follows from assuming that the atomic 
spins are all parallel to one another, which permits one to write Sj = (s/S)S 
for all i. The remaining sum equals J2i ^-ii which we call Itot, the total nuclear 
spin. For any given value of N n , I to t can take on values ranging in integer steps 
from N n /2 to either (if N n is even) or 1/2 (if N n is odd), with multiplicities 
that can easily be found. (See below.) Since Hm commutes with I 2 ot , we can 
consider the problem for one value of /tot at a time. Writing I instead of I to t 
henceforth, we arrive at a model Hamiltonian obtained by adding Eqs. (Q) and 
(§, i.e., 

H I = ^(-k 1 S 2 z +k 2 S 2 x -iu n S-l). (6) 

[For completeness, we give here the formula for the multiplicity, i.e., the 
number of times a given value of /tot appears when N n spins of magnitude 1 /2 
are added together. For / tot = (N n /2) — k, the multiplicity is given by 



k U-l 



(7) 



As an example, if N n — 6, multiplets with / to t = 3, 2, 1, and 0, occur 1, 5, 9, 
and 5 times respectively] 

It should be noted at this point that for a particle of diameter 50 A, which 
corresponds to S = 0(1O 4 ), and typical material or anisotropy parameters, the 
tunnel splitting Ao as given by Eq. (|J) is unobservably small. Of course the 
splitting goes up if 5* is decreased, but particles much smaller than 50 A seem 
difficult to attain controllably and reproducibly with present day technology.| 



b The last few years have seen some very interesting worl EjlOH on systems with much 
smaller values of S, of order 10. These systems are based on magnetic molecules, where the 
value of S is highly reproducible, but it is my belief that the interesting questions here are 
quite different, and have little to do with MQP. Interestingly, semicalssical methods, such as 
those developed in this paper, are likely to be quite valuable in analyzing these systems. For 
particles with intermediate values of S, say about 100, on the other hand, it will probably be 
necessary to devise new ways of treating the environment because the law of large numbers 
will no longer be applicable. 
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MQC is a more likely proposition in antiferromagnetic particle£jl2xE!J where 
the two states involved differ in the orientation of the Neel vector, or the^ciins 
on one of the sublatticesf] The essential aspects of the earlier papers, cm in 
which the influence of nuclear spins was studied using instantons, hold equally 
for ferro- and antiferro-magnetic particles. My goal in this article is to try and 
verify the predictions of the instanton approach as quantitatively as possible. 
It is simpler to do this using a ferromagnetic model, and I shall therefore limit 
myself to that in this article. 

I also do not wish to discuss at length the experimental observability of 
the reversal phenomenon. In addition to the values of the physical parameters 
mentioned above, this hinges on a number of other issues such as the tempera- 
ture, the signal size, stray magnetja-fields, the nature of the substrate, etc. We 
refer readers to previous papersffl€3E3 for detailed discussions of these points. 
For the purposes of this article, we will only note that the ratio Lo n /uj e is of 
order 10~ 3 to 1CP 2 in antiferomagnets, and is about one order of magnitude 
higher in ferromagnets. Throughout this article therefore, we shall assume that 
Lu n /uj e -C 1, and work to leading order in this ratio. We shall further assume 
that Ao/ujn <C 1, which, by virtue of the exponentially small WKB or Gamow 
factor, is almost certain to be the case for values of S of interest to us. 



2.2 What to Calculate; Preliminary Arguments 

Having settled on a theoretical model, let us ask what physical quantity we 
should calculate. One way the dynamical behaviour of the moment can be 
described is in terms of a time dependent probability tq_find it along some 
direction given an initial state in which it was preparedou This description 
best applies to experiments on a single system. Another way is to find the 
appropriate frequency-dependent dynamical susceptibility x( w )- This descrip- 
tion better applies to an assembly of identical or nearly identical systems, such 
as is obtained in an NMR experiment. We will adopt the second approach. 

Let us suppose then that we wish to calculate x"(u>) for our magnetic 
particle. Since the Hamiltonian (Q) describes a closed system, \" can on ly 
consist of a set of delta functions. The real system is of course not closed, 
and has means of energy relaxation such as phonons, with which the spins can 
interact in a variety of ways. These will lead to line widths and broadening 
in the usual way. Now in a setup such as NMR, the resonance frequencies 
are few in number and usually known quite accurately to begin with. Interest 

c The reason for this is quite simple. The tunnel splitting can still be written in the form 
The energy barrier Su is of the same order of magnitude, but the electronic spin attempt 
frequency, u> e , is about 100 times higher for antiferromagnets than ferromagnets, being given 
by the geometric mean of an exchange energy and an anisotropy energy in contrast to Eq. . 
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Figure 1: Structure of low lying energy levels of coupled electron and nuclear spin systems. 
The states | ± 0} are electronic spin ground states with S localized along ±z in the absence 
of the nuclear spins. The corresponding first two excited states are | ± 1} and | ± 2). The 
central part of the figure shows how the ground electronic levels are modified by the nuclear 
spins, still ignoring tunneling. The ratio u> e /oj n ^> 1 and is not accurately represented in the 
figure. If uj„ ^> Ao, only states joined by dashed lines mix with each other to any appreciable 
extent once tunneling is turned on. 



then attaches to the calculation of line shifts and widths and shapes and so 
on. In our problem, on the other hand, the presence of nuclear spins alters the 
resonance frequencies themselves drastically. To see this, let us denote by | ±0) 
the approximate eigenstates of the bare Hamiltonian Tio that correspond to the 
spin having a mean value (S z ) w ±S, i.e., to states in which the spin is localized 
in the ground state in one of the two energy minima centerd at ±z. Tunneling 
mixes these states and splits them by Ao. In addition to these states, let us 
consider the next few excited states in each well, and denote them by | ± 1), 
| ± 2), and so on. Since the small oscillations (or more accurately precessions) 
around the ±z directions are approximately harmonic, these states will form an 
approximately equally spaced ladder with spacing u) e , i.e., the states | ±n) will 
lie an energy oj e above | ± (n — 1)) for low values of n. (See Fig. 1.) Now let us 
bring in the nuclear spins. Since (±0|S| ± 0) = ±Sz to very good accuracy, the 
combined elctronic-nuclear spin states will be eigenstates of I z , which we label 
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by | + 0,p), and | — 0,p), where p — 2I Z , will be split into 21+1 Zecman levels 
with a spacing uj n . (The higher states | ± 1), | ± 2), . . ., will also be similarly 
split, but we do not consider that for the moment.) As long as / is small 
enough so that 2Iuj n <C w e , our basic assumption that u> n 3> Ao combined 
with the physically obvious but important fact that significant resonance is 
only possible between states that are degnerate to within the matrix element 
connecting them implies that the state | + 0,p) will resonantly tunnel only to 
the state | — 0, — p). Since this tunneling now involves a change in the nuclear 
spin state in addition to that of the electronic spin, the splitting should in 
general be different from Ao. We denote the magnitude of this splitting by 
A(I,p), i.e., 

A(I,p) = ±(E g(0tP) -E u{0:P) ), (8) 

where -E s (o, P ) aud -E«(o,p) denote the energies of the antisymmetric and sym- 
metric linear combinations (|0,p) =F \0,p))/^/2. Note that the antisymmetric 
state need not be the one with higher energy. 

To orient further discussion, let us summarize the results of the instanton 
approach. B'B In this approach too the starting point is that the nuclear spins 
spoil the degeneracy between the | ± 0) states. An instanton connecting degen- 
erate states must involve the flipping over of a certain number of nuclear spins 
along with that of S. Since the nuclear spins can only respond on a time scale 
uj~ 1 to any perturbation, and the instanton has a temporal width of order u~ , 
perturbation theory shows that each nuclear spin coflip reduces the tunneling 
amplitude.-hy a factor ~ (u; n /u; e ).g The amplitude for p nuclear spin coflips 
was foundcra (in the case of an antiferromagnetic particle) to be 

A p = (™ n /2w e )l p lA . (9) 

In other words, the tunneling amplitude decreases geometrically with the num- 
ber of units by which I z must change. If this is so (and we will find here that 
by and large it is), then clearly the more interesting problem is find the tun- 
neling frequencies A(I,p) themselves, and the associated spectral weight to 
be assigned to each frequency. The issue of linewidths and relaxation becomes 
secondary. As stated in the previous section, however, the calculation on which 
these conclusions are based is not totally satisfactory. It is with this viewpoint 
that we focus in this article on calculating the splittings as carefully as possible. 

d This reduction of tumbling amplitudes due to nuclear spin coflips was also found by 
Prokof 'cv and Stamp. til The highly chopped-up nature of the \" spectrum was missed 
by them, however, iWfd the mechanism for decoherence initially studied by them is rather 
different from mineLTu 



8 



It should be clarified that the amplitudes A p are not to be identified with 
the tunnel splittings themselves. To see this, consider a multiinstanton trajec- 
tory that starts and ends at a specific state with S z = S and I z = p/2. The 
intermediate states with a long residence time in this trajectory are those with 
S z = S, I z = p/2, and S z — — S, I z = —p/2. Since p spins can be chosen 
from N n in many ways, there can be several such intermediate states, and 
the amplitudes for all the corresponding trajectories must be added together. 
This leads to interference effects. When the combinatoric factors are added 
together, it is found that the splittings are given by El 



A p , n = \p\,\p\ + l,...,(N n + \p\)/2. (10) 



n 
>l, 

In this article we will try and see how correct these results are. The main 
flaw in the previous work is the conclusion that A p = A_ p , i.e., that the 
\0,p) <-> | — 0, —p) and |0, —p) <-* \ —0,p) splittings are the same. In yet more 
words, the splittings for p coflips are independent of whether the coflipping 
spins are parallel or antiparallel to the large electronic spin S. We will see that 
this is no longer true. Further, the geometric dependence (w„/(x» e )' p ' seems to 
be only approximately correct. Finally, we will show in Sec. 5 that the bino- 
mial factor in Eq. (|10| ) follows from a simple effective tunneling Hamiltonian 
connecting the states \0,p) and | — 0, — p), which holds as long as w n /w e < 1. 

One trivial point should be noted and disposed of once and for all. Since 
the Hamiltonian (|^) is time-reversal invariant, the bare tunnel splitting van- 
ishes unless S is an integer, and all the splittigns A(I,p) vanish unless S + I 
is an integer. For the purposes of interpreting the relation (|^) (or a more ac- 
curate replacement) it is useful to define Aq for all S by Eq. (|2|) or its more 
accurate version, Eq. (|24|) below. 

3 The Discrete WKB Method 

Like its continuous counterpart, the discrete WKBjHethod is applicable to a 
wide variety of settings. (See the review by BraunLL3 for examples.) We will 
only give an overview and physically motivated discussion of this method using 
the Hamiltonian (^) as an illustrative example. 

Let us write a general eigenfunction of Ho as J2 m a m\m), where as usual 
S z \m) = m\m). Schrodinger's equation then becoems a three-term recursion 
relation for the coefficients a m : 

W m a m + tm,m+2 a rn+2 + ^m,m-2 a m-2 = Ea m , (11) 

where w m = (m\T-lo\m), and t m ^ m ±2 = (TO|7io|w ± 2). Now for large S 1 , the 
differences w m +2 — w m and t„ ljm +2 — t m — 2,m, etc., are of order 1/S relative 
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m, m+ 2 




m — *~ 

Figure 2: Mapping of spin prohlfim onto an electron hopping on a lattice. The diagram is 
drawn for a case such as Eq. ( |ll| ) with only second neighbour hopping. It is evident that 
the sites connected by dashed and solid lines then belong to two disjoint subspaces of the 
Hamiltonian. 



to w m and t m;m+ 2 themselves. It is thus extremely useful to view Eq. ( |11| ) as 
arising from a tight-binding model for an electron on a one-dimensional lattice 
with sites labeled by to, on-site energies w m , and hopping energies t m rn ±2, 
that vary slowly with position.^ (See Fig. 2.) This viewpoint immediately 
suggests the approximation of semiclassical electron dynamics, and indeed this 
approximation is identical to discrete WKB. 

To apply semiclassical dynamics, we define local TO-dependent functions 
w(m) and t(m) by 

w(m) = w n 



(12) 

t(m) = (i m ,m+2 +tm,m-2)/2, 

which we extend to continuous values of m by demanding that they be smooth, 
and that dw/dm and dt/dm be of relative order S^ 1 with to formally re- 
garded as a quantity of order 1. A particle of energy E can be assigned a 
local, TO-dependent, wavevector q(m) in complete analogy with the continuous 
WKB approach. The only difference is that the kinetic energy is given by 
2t(m) cos q(m) instead of <7 2 (m)/2/Lt (/z being the mass). We thus obtain 

/ \ i / E — w(m) \ 

«= cos (r^r)' (13) 



The general solution to (11) is given by linear combinations of 



1 ( f n , ,,dm'\ 

:exp ±t / q(m')—\. (14) 



e For the Hamiltonian ill, the hopping connects sites differeing by Am = 2, so that the 
eigenvalue problem divides into two subspaces depending on whether S — m is even or odd. 
This point is of no consequence for application of the discrete WKB method itself, as the 
problem can be recast as one of nearest-neighbor hopping in each subspace. 
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Figure 3: Discrete WKB method potential energy diagram for a symmetric double- well 
problem. The functions U±(m) = w(m) ± 2t(m) are the local band edges. For a particle 
with energy E as shown, the regions mi < m < m,2, and 1T13 < m < mn are classically 
accessible, and represent the two wells. The particle must tunnel across the classically 
forbidden region between 7713 and m^. 



In this equation, 

v(m) = — 2t(m) sinq(m) (15) 

is the local particle velocity (equal to dE/dq), and the factor kv(m)] -1 / 2 fixes 
the normalization so that the probability current is conserved. Also, the factor 
of 1 /2 in the integrand reflects the step length Am = 2 in Eq. ( pi] ) . 

We conclude here our general discussion of the discrete WKB method, 
and continue the illustration with the Hamiltonian (^) in the next section. 
Braun's reviewli3 contains a proper proof of the above aprroximations, as well 
as connection formulas at turning points, Bohr-Sommerfeld quantization rules, 
etc. Most of these are physically apparent, and the chief novelty arises from 
the fact that for particle in a periodic potential (or a tight-binding model), the 
allowed energies lie in a band, and are bounded both above and below. This 
gives rise to turning points when the energy equals the local upper band edge 
w(m) + 2\t(m)\i in addition to those that occur when the energy equals w(m) — 
2\t(m)\. These features are illustrated in Fig. 3. Braun refers to the first kind 
of turning point as "unusual" , but the associated formulas are easily derived 
from those for a "usual" turning point by making a gauge transformation which 
changes the sign of every other coefficient a m . The associated bookkeeping can 
be fairly cumbersome, however, and it is easy to make mistakes. 
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4 Tunnel Splitting for the Hamiltonian ([!]) 

We now turn to applying the discrete WKB method to finding the ground 
state wavefunction and tunnel splitting Ao for the_Hamiltonian ([[]). (We can 
also find the higher spittings and wavefunctions. cj) There are three steps in 
finding the wavefunction itself: (i) find it in the classically allowed region, (ii) 
find it in the classically forbidden region |m| <§; S, (iii) match the two parts of 
the wavefunction using the connection formulas or otherwise. The splitting is 
then obtainable by a textbook formula.^] 

To find the wavefunction in the classically allowed region near m = — S, 
e.g., it is advantageous to write C n = a_s+„, with n = 0, 1, . . .. It is generally 
necessary to obtain the functions w(m) and t(m) to an accuracy such that the 
first two terms in an expansion in powers of 5 _1 are correctly given; if only 
the leading term is kept, then retention of the (v(m)) factor in Eq. (14) 



can not be justified. In carrying out this exercise near the ends of the chain, as 
is the case now, it is necessary to treat n = S + m rather than m as a quantity 
of order S°. It is also useful to add a constant fc^S'fS' + 1) to the Hamiltonian. 
The Schrodinger equation then reads 

(2fci+fc 2 ) (n+ C n + ht 2 ^n(n - l)C„- 2 + ^ W(™ + l)(n + 2)C n+2 = EC n 

(16) 

This equation, however, can be solved exactly. It is that for the harmonic 
oscillator Hamiltonian 

H ho = (2fci + k 2 )a)a + h 2 (a 2 + a^), (17) 

where and a are raising and lowering operators for the number eigenstates 
\n) obeying a) a\n) — n\n), and C n — (n\ip) for an energy eigenfunction \ip). 
We can diagonalize Hho by writing a — (x + ip)/2 1 / 2 , = (x — ip)/2 1 / 2 , where 
x and p are canonical position and momentum operators obeying [x,p] = i as 
usual. This yields 

Hho - k 12 x 2 + k lP 2 . (18) 



* See Eq. (03) below. This formula is anologous to one thaf-iappears in the solution to 
Problem 3, sec. 50, in the famous text by Landau and Lifshitzcj However, it is not the final 
formula in that solution! We have found it preferable to use the unnumbered intermediate 
equation in which the splitting is given as a product of the wavefunction and its derivative 
at the symmetry point of the potential. The three dimensional analog of the latter formula 
is often named after C. Herring who first used it find the splitting of the electron terms in 
the H^" ion. This problem is also discussed by Landau and Lifshitz in the solution to the 
problem accompanying Sec. 8f . 
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The wavefunction C n (automatically normalized to unity) is now easily ob- 
tained by evaluating the overlap {n\if)) in the x representation. We quote the 
final result for the ground state. C n vanishes for odd n, and for even n it 
equals 

1/2 / tanh^ n/2 



c n = V^re^ n/2 ) , (19) 

where 

6 = (l/4)ln(fei a /fei). (20) 

It is also useful to note the relations sinh2# = cosh 28 — (2k\ + k 2 )/u> e , 

and tanh26> = k 2 /(2k 1 + k 2 ) 1 with u> e = 2{k 1 k 12 ) 1/2 as defined in Eq. (|). That 
uj e is the small precession frequency for the electronic moments is now obvious 



from Eq. (|1S|) 



[We mention for completeness that the above results also follow from mak- 
ing a Bogoliubov transformation 

a = cosh0 b — sinh0 w , . . 

at = -sinh0 6 + cosh (9 &t, ( 21 > 

which, with 6 given by Eq. (|20|), diagonalizes 7iho- The coefficients C n follow 
from an obvious iwo-term recursion relation connecting C n and C n _2-] 

The only constraint used in deriving Eq. ( |l9| ) was n <C S. In particular, 
the result holds for S 3> n ^> 1. Thus it actually extends out into part of the 
classically forbidden region. This is a great advantage as it allows us to match 
to the quasiclassical form of the wavefunction on the other side of the turning 
point directly, obviating the need for connection formulas. 

The wavefunction in the classically forbidden region is directly found using 
the formulas (|l4| ) and (|l5|). Under the conditions n/2 ^> 1, S — n/2 1, we 
find that for even n 

C n = (tt| (l - ^) J (sech^-tanhfl)"/ 2 . (22) 

(C n continues to vanish for odd n.) We have already adjusted the multiplicative 
constant in this result, so as to agree with Eq. (|l9|) when S 3> n^> 1. 

The tunnel splitting is generally given by the formula (see footnote /) 



2t 2 ,o\a Q {a 2 - o_ 2 )|, even S, 
2i 1 ,_ 1 |af — (a_i) 2 |, odds'. 



(23) 



(It is also easily shown that the antisymmetric state is lower in energy for odd 
S, and higher for even S.) Applying this formula to the problem at hand, we 
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Table 1: Comparison between numerical and analytical [Eq. (124)] results for the tunnel 
splitting without nuclear spins. The parameters are k\ = 5.0, = 20.0. 



s 


Aq (numerical) 


Aq (analytic) 


Error(%) 


10 


9.282 x 10" 


3 


9.749 x 10" 


3 


5.0 


11 


3.738 x 10- 


-3 


3.906 x 10" 


3 


4.5 


12 


1.497 x 10" 


-3 


1.558 x 10" 


3 


4.1 


13 


5.974 x 10- 


-4 


6.195 x 10" 


4 


3.7 


14 


2.375 x 10" 


-4 


2.455 x 10" 


4 


3.4 


15 


9.412 x 10" 


-5 


9.708 x 10" 


5 


3.1 


16 


3.721 x 10" 


-5 


3.830 x 10" 


5 


2.9 


17 


1.468 x 10" 


-5 


1.508 x 10" 


5 


2.7 


18 


5.778 x 10" 


-6 


5.927 x 10" 


6 


2.6 


19 


2.271 x 10" 


-6 


2.326 x 10" 


6 


2.4 


20 


8.910 x 10" 


-7 


9.115 x 10" 


7 


2.3 



obtain 



S 



1/2 



A = 4w e ( - I sech 6 tanh s 6. (24) 



Since tanhf? < 1, the last factor in this result can be rewritten as a WKB 
or Gamow factor exp(— S\ lntanh#|). For k\/ki2 <C 1, this is approximately 
equal to exp(— ASki2/u> e ), which is of the same form as Eq. (|J). 

It is of interest to note that the only previous calculations of Ao to the-same 
accuracy (and with which we agree completely are by Enz and Schilling, lj and 
by Belinicher, Provdencia, and Providencia,E3 both of whom employ methods 
of much greater corac 



dty Other tunnel splitting calculations either find an 
incorrect prefactor, Ej'Ea or not at all. El In Table 1, we compare the formula 
with exact answers for the splitting obtained by numerical diagonalization of 
Ho for S ranging from 10 to 20, for parameters &2 = 20.0, and ki = &2/4 = 5.0. 
This choice yields uj e = 22.36, tanhtf = 0.3820, and sechf? = 0.9242. As can be 
seen, the errors are quite small and the results improve steadily with increasing 
S. Note further that Eq. ( p4] ) consistently overestimates-the splitting. These 
results are in keeping with previous numerical studies. EjeZ 

5 Inclusion of Nuclear Spins 

5.1 One- Co flip Splittings 

The results of the previous section enable us to calculate the tunnel splitting 
A(J,p/2) for small values of I and p. Let us begin with the simplest possible 
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777 



Figure 4: Lattice diagram for MQC problem with one nuclear spin cofiip. The k x S^. and 
S±I+ bonds are shown by heavy solid and dashed lines respectively. The division into two 
subspaces is again eident. 

case, when / = 1/2, and p = ±1. Then as discussed in Sec. 2, it is necessary 
to take S to be a half integer. We will write S = N + 1/2, where N is an 
integer. 

In the language developed in the last section, the problem is now equiva- 
lent to that of an electron hopping on a lattice with two rows, corresponding 
to I z = ±1/2, and 2S + 1 = 2(N + 1) lattice points per row. The problem is 
represented pictorially in Fig. 4, where the various non-zero off-diagonal ele- 
ments of the Hamiltonian (^|) are shown as bonds connecting the corresponding 
lattice points. The bonds connecting alternate sites in the same row arise from 
the k2S%. term, while the bonds connecting sites in different rows arise from 
the terms S+I— and S—I+. It is apparent that the sites form two groups that 
are disconnected from each other, or in proper mathematical language, that 
the Hamiltonian divides into two disjoint subspaces. In particular the split- 
ting of the two lowest states in the space that contains the state^] | —5,-1/2) 
equals A(l/2, 1), while the splitting in the space containing | — S, +1/2) equals 
A(l/2,-l). 

Let us first calculate A(l/2, — 1) [which equals A_i as we will show in 
Eq. j46|) below]. The lattice diagram for this is shown in Fig. 5, where we have 
thrown away the states in the other subspace, and we have relabeled the sites 
as explained below. We will obtain the splitting by finding a suitable effective 
two-state Hamiltonian. To this end, let us first ignore the S+i_ and S-I+ 
terms in Hi- The two lowest energy states are then | + 0,p = — 1) = |I) and 
| — 0,p = +1) = | II). These states live on different rows of the lattice, with 
1 1) and | II) being localized near the right edge of the lower row, and the left 

9 It is necessary in this section to refer to two different types of states with similar labels: 
the states \S Z ,I Z ), and the states | ±n,p), where n is the harmonic-oscillator-like excitation 
index used in Sec. 2.2. We shall avoid notational confusion by refering to the second type of 
state as |0, p = 1), | — 0, p = 1), etc., in which the label p is explicitly identified. 
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1 2 J N-l N 

* * • // • * 1 * 

\ \ \ \ \ 
\ \ \ \ \ 
• • • • // » • • 

N N-l 2 10 



Figure 5: Lattice diagram and labelling scheme for calculation of overlap for A(l/2, — 1). 
N = S — 1/2. Only sites in the subspace including Sz = —S, Iz = 1/2 are shown. 



edge of the upper row, respectively. It should be carefully noted that neither 
state is large near both ends. The state which is large near S z = —5+1 
and I z = —1/2 (lower left corner), for example, is approximately equal to 
| — l,p = —1) which is higher in energy than |I) or |II) by u) e . The additional 
bias in the on-site energy due to the term u n S z I z /S is not large enough to offset 
this energy difference since w„ < w e . For the same reason, the wavefunctions 
of the states |I) and |II) are very well approximated by the wavefunction found 
in the previous section. It is necessary only to redefine the position coordinate 
appropriately. We do this by labelling the sites in Fig. 5 by j = 0, 1, . . . , N, 
where j runs from left to right for the upper row, and in the opposite direction 
for the lower row. Then 

(S-2j, -1/2|I) = (S + 2j,+l/2\ir) = C' j . (25) 

Further, Cj = C2j where Cij is given by Eqs. (|l|) and @. 
The effective two-state Hamiltonian can be written as 

n eS =(^ I), (26) 



V E 

where E a = (I\Hi\T) = (II|Wj|II), and V = (I|W/|II) = (H|H/|I). The value of 
Eq is clearly immaterial, while the tunnel splitting is given by 2\V\. The off- 
diagonal element V arises from the S+I- and S-I+ terms and can be directly 
evaluated because w„ < w e . We have 



N I 1 

V = Y,C' j C' N _J-S + 2j + l,— 



2S S+I - 



J=0 

JV 

= -K/2S) C'f' N _, } [{2j + 1)(2S - 2j)] 1/2 . (27) 

3=0 
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N N-l J 2 1 

* * • // • * * • 
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\ \ \ \ 
• • * %-H- — • • • 
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Figure 6: Same as Fig. 5 for A(l/2, 1). Only sites in the subspace Sz = S, I z = 1/2 are 
shown this time. 



The dominant contribution to the sum arises from the middle of the chains, 
i.e., from values of j such that (N — ^> 1. Using Eq. ( p2|) and neglecting 
terms of order unity in comparison with j and N — j, we obtain 



V = -(-l) N 2uj n {S/Tr) 1 / 2 sccheta,nh N eJ^iU + l/2)(# - j + 1/2)] 



-1/2 



(28) 

For JV> 1 the sum may be approximated by an integral, which is easily shown 
to equal tt. Making use of Eq. (pi|), and recalling that S = N + 1/2, we can 
write the tunnel splitting (which equals 2\V\) as 

A_i = A(l/2,-l) = - — (tanh6>)- 1/2 A . (29) 

2 LU e 

The splitting A(l/2, 1) (= Ai) can be found in the same way. The ap- 
propriate lattice diagram is drawn in Fig. 6. This time the two low energy 
states are \ + 0,p = 1) and | — 0,p — —1), and the wavefunctions are given 
by the same C'j as above with the sites labeled as shown. The key quantity is 
again the off-diagonal element. It is obvious from the figure that this is again 
given by a sum as in Eq. (|27|) . The only change required is that we replace the 
product of wave functions in that summand by C'^C' N+l .. The sum is again 
dominated by the central region, and since C*j +1 ~ — tanh in this region, 
the change in the summand has the effect of multiplying the each term in the 
sum by — tanh#. The final result for the splitting is, therefore, 

Ax = A(l/2, 1) = ^(tanh^^Ao. (30) 

2 co* 



We thus see that apart from the factors of (tanh6') ±1 ^ 2 , the results (29) 



and (|30D agree with the instanton answer, Eq. (H). This agreement shows that 
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Table 2: Comparison between numerical and analytical [Eqs. ( pj| ) and (pc|)] results for onc- 
coflip tunnel splitting. The parameters are k\ = 5.0, ki = 2CT0, and io n = 0.1. For each 
value of S, the entry in the upper row is Ai, and that in the lower row A_i. 



s 


A±i (numerical) 


A-j-i (analytic) 


Error(%l 


10± 

2 


2 576 x 10" 


5 


2.681 x 10" 5 


4.1 




6.661 x 10" 


5 


7.018 x 10~ 5 


5.4 


2 


1 034 x 10" 


5 


1.072 x 10" 5 


3.6 




2.674 x 10" 


5 


2.805 x 10~ 5 


4.9 


2 


4.134 x 10" 


6 


4.267 x 10~ B 


3.2 




1.069 x 10- 


5 


1.117 x 10" 5 


4.5 


13± 


1.646 x 10- 


6 


1.694 x 10" e 


2.9 




4.257 x 10" 


6 


4.435 x 10~ 6 


4.2 


14± 


6.534 x 10" 


7 


6.705 x 10-'' 


2.6 




1.689 x 10" 


6 


1.755 x 10~ 6 


3.9 


15± 


2.587 x 10" 


7 


2.648 x 10" ' 


2.4 




6.688 x 10" 


7 


6.933 x 10" 7 


3.7 


16± 


1.022 x 10- 


7 


1.044 x 10~ 7 


2.2 




2.641 x 10- 


7 


2.732 x 10" 7 


3.5 


17± 

x 1 2 


4.024 x 10- 


8 


4.105 x 10~ 8 


2.0 




1.041 x 10" 


7 


1.075 x 10" 7 


3.3 


18± 


1.583 x 10" 


8 


1.612 x 10" 8 


1.9 




4.092 x 10" 


8 


4.221 x 10~ 8 


3.2 


19± 


6.217 x 10" 


9 


6.322 x 10" y 


1.7 




1.608 x 10" 


8 


1.655 x 10" 8 


3.0 




2.438 x 10" 


9 


2.476 x 10~ 9 


1.6 




6.303 x 10" 


9 


6.482 x 10" 9 


2.9 



the instanton method is basically sound, and although it can not in its sim- 
plest version give the prefactors in Ao correctly, the overall picture it provides 
is correct. In fact, the (tanhf^ 1 / 2 factors can also be found from the instan- 
ton approach by exercising a little more care. These origin of these factors 
is actually very easy to understand. They reflect the fact that the state with 
antiparallel alignment of the electronic and nuclear spins has greater quantum 
fluctuations or zero point motion than the state with parallel alignment. An- 
other way to say this is that the state \S Z = S, I z = I) is an eigenstate of the 
hyperfme coupling term in Tix, while \S Z = S, I z = — I) is not. This difference 
in zero point fluctuations shows up as a higher tunneling rate for coflips in 
which the nuclear spins are oppositely aligned to the electronic ones. 
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12 n-1 n n+1 

Figure 7: Lattice diagram for the case 1=1. The bonds entering into the calculation of the 
transition matrix element V relevant to A(l, 2) and A(l, —2) are shown as dashed and solid 
diagonal lines respectively. 



We also remind readers that in Eqs. fl29[ ) and Q3CQ , Ao is not the tunnel 
splitting for the bare problem without nuclear spins (which vanishes because 
S is half- integral), but rather the pseudo-splitting formally given by Eq. (p4|). 

In Table 2, we show the comparison between the above formulas and exact 
answers from numerical diagonalization of Hi with I = 1/2. We choose ki — 
5.0 and ki = 20.0 as before, and uj n = 0.1. Since u e — 22.36, the condition 
Lj n <C w e is well obeyed. Again the agreement is very satisfactory and improves 
with increasing S. In particular, the inclusion of the (tanh#) ±:L ' 2 factors seems 
to be required. 



5.2 Two-Coflip Splittings 

Let us proceed to calculate the two-coflip splittings, A(l,±2) (= A±2), by 
studying the Hamiltonian (0) for I = 1 and integer S. As before, we carry out 
the calculation by including the S Z I Z term in the unperturbed Hamiltonian 
along with Ho, and treating the S+I— and S—I+ terms as a perturbation. 

Suppose we wish to find A(l, 2). The appropriate lattice diagram is shown 
in Fig. 7. This tunneling process requires making a transition from the state 
| — 0,p = —2) = |I) to the state |0,p = 2) = |IV). However, there is no 
direct matrix element of Hi between these states, and we must go through 
intermediate states. The lowest energy intermediate states are | — l,p = 0) = 
| II) and \l,p = 0) = |III). Keeping only these we obtain a four-state effective 
Hamiltonian 

Weff = n Adlrt t? ' t/ I ' ( 31 ) 
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where the equality of various elements is assured by symmetry. Here, Eq = 
<I|W/|I) = (TV\Hi\TV), Ei = (ll\Hi\ll) = (III|Hj|III}, V = (H|Wj|I), and 
A' 1 ' is the tunnel splitting between the states |II) and | III) - This splitting is 
identical to the splitting between the firsi_fixcited states of the Hamiltonian 
Tio without nuclear spins, and is given bycj 

A<,) - ;d£sh? a »' (32 > 

Secondly, the energy difference Ei — E = ui e + 0(u) n ). 

The key quantity in calculating A2 is thus the matrix element V as before. 
This is given by a sum involving the lattice site wavefunctions of the states 
1 1) and I II) and a J_5+ matrix element. Refering to Fig. 7, consider the 
contribution of the bond connecting site n in the lower row to site n — 1 in the 
middle row. The matrix element for this bond equals 

-^(-S + n-lJ z = 0\S + I-\-S + nJ z = -l) = -^(n(2S-n + l)) 1/2 . 

(33) 

The wavefunctions at the two ends of this bond, on the other hand, are given 
by (n\ipo) and (n — l\ipi), where | Vo) an d |^i) are the ground and first excited 
state of the harmonic oscillator Hamiltonian Tiho, Eq- (p"7|), and \n) are the 
eigenstates of the number operator a^a. Thus, 



V = 



-g- J2H 2S ~ n + !)) <^|n - l)(n\ih) 

OO 

v n— 

00 



= ^=smh6. (34) 

The second equality in this chain follows from noting that the sum is dominated 
by values of n of order unity, which allows us to neglect n in comparison to S. 
The final equality follows if one makes use of the Bogoliubov transformation 
(|U). Then, (V>i|a|V>o) = - sinh6»(^i \tf\ip) =-sinh0. 
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Table 3: Comparison between numerical and analytical results for two-cofiip tunnel split- 
tings, for k± = 5.0, &2 = 20.0, u> n = 0.5, S = 15, and 1 = 1. In the column labeled 'ratio' we 
give the ratio of the numerically obtained splittings to those given by the instanton method, 
i.e., to Eqs. Q and (Ei), but with a numerical factor of 7r 2 /4 instead of 2. 



Numerical Ratio 
2.559 x KT 8 0.56 
A_2 3.488 x 1(T 7 1.11 



To leading order in u> n , the energy splitting of the lowest two states of the 
effective Hamiltonian ([H]) is easily shown to be given by 

V 2 V 2 V 2 m 

7 A (1 \ (35) 



E l -E - A« /2 E x - E + A« /2 (Ex - E f 

Combining Eqs. (^), (p^), and using the fact that Ex — E Q = uj e + O(oj n ), we 
finally obtain 

LU 2 

A 2 = A(l,2)w 2-£ tank 0A O . (36) 

The calculation of A(l, —2) proceeds very similarly. The effective Hamil- 
tonian has the same structure as Eq. ([H]), except that this time the states |I) 
and | IV) must be taken as | — Q,p = 2) and \ +0,p = —2) respectively, while 
the states |III) and |IV) are unchanged. This time the bonds connect site n in 
the upper row with site n + 1 in the lower row, and the matrix element is 

- ^(S + n + 1,1, = 0\S_I + \ - S + n,I z = 1) « -^=(n + l) 1 / 2 (37) 



instead of Eq. (]33|). Using the same technique as before, the transition matrix 
element V is found to be 



V = ^2(n + l) 1/2 (ipx\n + l)<»|Vb) = -S- 1/2 uj n coshf?. (38) 



^ n=0 

Substituting this in Eq. (|35|), we obtain the splitting as 



U] 2 

A_2 = A(l,2) Ps2-£(tanh0)- 1 A o . (39) 



The results ( |36| ) and (39) are very close to but not exactly what one would 
expect from the instanton method, even after the fluctuational factors of tanh 9 
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have been included. If we define rj± — (7ru; n /2a;e)(tanh0) : ' =1 / 2 , then an instan- 
ton argument would lead us to expect 

A± 2 = 77±A . (40) 

This would require the numer ical factors in Eqs. @ and (||) to be tt 2 /4 
instead of 2. Could this really be the case, and if so, what is the source of the 
discrepancy in our present calculation? One possible answer is that we have 
ignored the contribution of higher energy intermediate states in the calculation 
of A±2- In fact, it can be shown that these states' contribution is also formally 
of order (cj„/cj e ) 2 Ao, but with different numerical prefactors, which decrease 
rapidly with increasing intermediate state energy. The accurate calculation of 
these numerical factors is difficult. Similarly, we neglected the terms of higher 
order in S^ 1 in expanding the matrix elelment (f33|). Terms of such higher 
order are also present in the contributions from the higher intermediate states, 
and it is not clear that the sum of all these terms will continue to be of higher 
order in S . When these approximations are taken into account, we cannot 
exclude the possibility that Eq. ( |40| ) holds exactly. 

We show a comparison between our theoretical and numerical results in 
Table 3. The quality of agreement is still fairly good, though not as impressive 
as for the one-coflip splittings. We have not carried out this numerical work as 
extensively as in the previous cases, and so can not comment on the behavior 
with increasing S. 

5.3 Higher Co flip Processes; Effective Hamiltonian 

It is evident that the discrete WKB method is increasingly ill suited to the 
calculation of higher coflip processes. Systematic numerical investigation of 
this problem runs into the difficulty that the splittings necessarily decrease 
with increasing coflip number p, requiring the numerical diagonalization to be 
carried out to increasing precision. My own numerical studies are very limited, 
and while it is undoubtedly possible to improve on them, that would require 
substantially greater time and effort than I am prepared to commit! 

There is, however, one aspect of the higher coflip problem, which can be 
understood quite simply. [We have alluded to this earlier; see the discussion 
immediately preceding Eq. ([To|).l This is that splittings A(I,p) for different 
/ but equal p can be related to each other. To see this, let us consider the 
problem in terms of N n nuclear spins each with spin 1 /2 rather than the model 
(||). There are then a large number of states in the | + 0,p) group that are 
degenerate with one another, and with an equally large number of mutually 
degenerate states in the | — 0, —p) group. (Recall that p = 2/' ot .) A p-coflip 
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process can take us from any of the states in the first group to several states in 
the second group. A convenient algebraic method for keeping track of which 
states are connected is as follows. To avoid cluttering up the formulas, it is 
best to treat the cases p > and p < separately. Let us do the p > case 
first. Let ir± p be projection operators onto the set of states with J* ot = ±p/2, 
irrespective of the value of I. Further, let 

<r+ = |+0>(-0|, 

a_ = |-0>(+0|, (41) 



be transition operators for the large electronic spin between the two states 
involved in the MQC process. Finally, let A p be the amplitude for the transi- 
tion between any two of the states in the groups with opposite I z . Then, an 
effective Hamiltonian that describes all the coflip processes of order p is 

Hf = A p {a + Tr p Q+K- p + h.c.) ; (42) 
h, J2,...,j p 

where = I x ± il y as usual, and the sum in Eq. ( |43] ) is over all distinct 
p-tuplets of indices chosen from the indices 1,2,..., N n . 

Let us check that Eq. ( (42"| ) does what it is supposed to. Consider the 
first term on the right, and let it operate on a ket. The projection operator 
7T_ p ensures that this term is only relevant when we operate on a state with 
I z = —p/2. The operator Q+ then raises (or flips) p nuclear spins leading to 
a state with I z = p/2. The projection operator n p in this term is added to 
ensure that the result is sensible when we operate on a bra. The second term 
in (42) ensures hermiticity and describes processes in which I z is lowered. 

It is now very easy, however, to find the / dependence of the matrix ele- 
ments of Hp { . We first note that because {I^) 2 = for any j, we can write 

Qi = ^ {it + it + ■ ■ ■ + ik) P = ^ + r- (44) 

The only non-zero matrix element of Tt p f between states of given / is thus 
equal to 



A„ /. 1 



i \ a ri+y 



p f(l,m=-p(I + yi,m=--pj=A p ^ (45) 

The case p < is now easily treated by minor changes of notation. We 
simply change Qp m the first term in Eq. Ji^ ) to its hermitean adjoint Q~ , 
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Table 4: Numerical results for higher coflip tunnel splittings. The parameters are fci = 5.0, 
k 2 = 20.0, w n = 0.5, S = l^and / = 0,1,2,3. The table shows A(J,p) divided by the 
combinatoric factor in Eq. (kffl. This combinatoric number is shown next to the result in 
parentheses. 



p 


/ = 


I = 1 




1 = 2 




/ = 3 







9.412 x 10" b (l) 


9.435 x 10- 


b (l) 


9.479 x 10" 




9.546 x 10" 


~ b (l) 


-2 




3.488 x 10- 


7 (1) 


3.497 x 10" 


- 7 (3) 


3.507 x 10" 


" 7 (6) 


2 




2.559 x 10" 


«(!) 


2.586 x 10" 


- 8 (1) 


2.628 x 10" 


-8(6) 


-4 








2.134 x 10" 


-9 (1) 


2.138 x 10" 


" 9 (5) 


4 








2.655 x 1Q- 


"(1) 


2.580 x 10- 


U (5) 


-6 












1.735 x 10- 


"(1) 


6 












3.624 x 10- 


13 (1) 



with the sum in the analog of Eq. (|43|) being over all |p|-tuplets. Thus Q p = 
(7~)|p|/|p|! for p < 0, and the matrix element analogous to ( p5| ) is identically 
evaluated. The upshot is that we can write 

for all p, positive or negative. 

We show in Tables 4 and 5 some numerical results for tunnel splittings 
involving up to 6 nuclear spin coflips. Instead of tabulating A(I,p) itself, we 
have divided it by the combinatorial factor in Eq. (|46"|). We expect that the 
resultant quantities will be independent of I, and we can see from the tables 
that indeed they are. In all the cases, the spread in the values of A p so deduced 
is less than 2%. 

Table 5: Same as Table 4 for odd numbers of coflips. Now S = 31/2, and 7 = 1/2, 3/2, 5/2. 
The other parameters are unchanged. 



p J=1 / 2 . 1 = 3 / 2 .. 1 = 5 / 2 . 

7L 3.265 x 10- y (l) 3.272 x 10^(2) 3.284 x 10~ e (3) 

1 1.327 x 10" 6 (1) 1.326 x 10" 6 (2) 1.323 x 10~ 6 (3) 

-3 1.619 x 10~ 8 (1) 1.623 x 10~ 8 (4) 

3 2.305 x 10~ 9 (1) 2.306 x 10~ 9 (4) 

-5 1.136 x 10~ 10 (1) 

5 6.457 x 1Q- 12 (1) 
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The fact that A(I,p) grows with I may make one wonder if this might not 
be a way to boost the resonance frequency for some of the tunneling processes. 
In fact, the result (f|6|) only applies to the ideal case where the hyperfine cou- 
plings are identical for all the nuclei. This limits the relevance of this result 
for actual systems to small values of p or /. Real particles typically possess 
some small spread in these couplings. This spread spoils the degeneracy of all 
the states in the \0,p) group (and likewise for the states | — 0, —p)). When the 
spread starts to equal A p , then the constructive interferences which give rise to 
the large combinatoric factors in Eq. ( flif ) are no longer possible. The spectral 
weight in x" is then shifted from the higher frequency peaks in a given cofiip 
group to lower frequencies. Secondly, we have neglected incoherent processes 
involving the nuclear spins. MQC is destroyed by a single such process, and 
since the likelihood that at leat one such process will occur in a given time 
grows linearly with N n , the number of nuclear spins, it is evident that it does 
not pay to increase this number. 



6 Conclusions 

The main purpose of this article has been to verify previous calculations Bi 
of the cofiip tunneling frequencies without using instanton methods. We find 
that by and large, the instanton calculations lead to correct answers. 

Our discussion up to this point has treated the system as closed, with no 
means for energy to flow in and out. If this were strictly true, this would imply 
that x" consisted of a sum of delta functions at the frequencies A(I,p). This is 
of course an idealization. To complete the discussion of the problem, we must 
also consider relaxation. If we put in broadening by hand, the qualitative pic- 
ture of the x" spectrum is as shown in Fig. 8. A proper investigation of the 
physical mechanisms behind this would take us into very different territory. 
Two obvious mechanisms which may be mentioned in passing are phonons 
and dipolar magnetic helds due to other nuclei such as protons which may 
be present in the material. It does not seem worthwhile to develop a quanti- 
tative theory of the relaxation due to even these processes in the absence of 
experimental impetus, but certain broadly valid comments can nevertheless be 
made. By appealing to the general principles laid down in Refs. 11 and 29, it 
is quite plausible that the relaxation can be treated by coupling the resonating 
variable (which we would regard as different for each cofiip line) to a phe- 
nomenological harmonic oscillator bath (which would also be different for each 
cofiip line). Broadly speaking, such a bath has two effects. First, it leads to a 
pulling down or downward renormalization of the bare coherence frequency A& 
to Aft. At the simplest level this renormalization can be understood in terms 
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Ipl = 2 



p = 



p =i p^- 1 



log CO 



Figure 8: Sketch of the expected tunneling susceptibility spectrum \" (u>) in the presence 
of nuclear spins, showing the first few polarization blocks. Note that the frequency scale 
is logarithmic. We have tried to indicate the fine structure within the higher polarization 
blocks, but we have made no attempt to represent the broadening of the lines correctly. 

of a multiplicative Franck-Condon factor. Second, it gives rise to true damp- 
ing of the resonance, i.e., an intrinsic broadening or linewidth. (There may of 
course be additional inhomogencous broadening arising from a spread in Af,.) 
The general point which is noteworthy is that as the bare frequency A& of 
the coherence phenomenon under study decreases, the environment suppresses 
it even further. This is especially so if the environment has any subohmic 
or ohmic component. Environmental degrees of freedom that might not have 
been relevant if Ah had been higher by a factor of 10, say, gang up on the 
system, as it were, and do start to matter. Further, degeneracy is broken by 
weaker and weaker stray fields and drifts, effectively eliminating the resonance 
altogether. 

All this implies that the higher coflip resonances in our problem are highly 
likely to be overdamped and unobservable. This point can be seen even in 
the simplistic model introduced in Rcf. 9, and briefly touched upon at the 
end of the previous section. If a single nuclear spin undergoes an incoherent 
transition from I z = 1/2 to 1 ~ z = —1/2 or vice versa at a rate 1/t, then any 
MQC in the particle as a whole is damped at a rate N n /r. All the lines in 
thus acquire a width of order N n /r irrespective of the number of coflips 
involved, and the higher coflip lines with A(I,p) < N n /r essentially merge into 
a mushy continuous background with no clear feature that could be identified 
as a resonance. This in turn means that while an accurate calculation of 
A(I,p) for large \p\ would still have theoretical interest, the unlikelihood of its 
relevance to experiment reduces the urgency of this calculation considerably. 
The most interesting remaining problem at this stage appears to be that of 
broadening of the very low order coflip lines. 

Let us therefore summarize our conclusions about how MQC in small mag- 



2G 



netic particles is affected by nuclear spins. The effect is severe, and the spec- 
trum of x" is broken up into a large number of lines which can be grouped by 
the number of nuclear spin coflips required to maintain degeneracy, or more 
conveniently by the nuclear spin polarization 



where Sj and Ii are the electronic and nuclear spins on atom i. The frequency 
of a line decreases quasi-geometrically with the number of coflips \p\, and there 
is a further distinction between p > and p < 0, with the latter having a 
higher frequency. Within each polarization group, there is a fine structure to 
the tunneling spectrum controlled by interference between the different ways 
in which p nuclear spins can coflip. The total spectral weight in a given polar- 
ization block is therefore given by f p , the probability that a particle will have 
polarization p. This probability can be controlled by thermostatstical factors. 
For example, if we assume that the nuclear spins arc in equilibrium with the 
elctronic spin, then the Boltzmann weights for 1^ parallel and antiparallcl to 
Si are in the ratio 1 : exp(— (3u) n ), where (3 — 1/ksT. This leads to 



where p = N n tanh(/3w n /2), and a p — N n sech (/3w n /2). The dominant reso- 
nance line is thus that with zero polarization, and has a weight fo- 

Finally, we note that the effect of nuclear spins on magnetic particle MQC 
is completely different from that of an oscillator bath. ISeI While unlike the 
latter, the present problem does not seem to lend itself to the study of elegant 
mathematical questions such as renormalization group connections with the 
Kondo problem, t2rtil it provides a novel and interestingly different example of 
how interaction with the environment suppresess MQC. 
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(47) 



f p = (2™l)-V*e-b>-W< 



(48) 
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